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A random-phase approximation (RPA) treatment of edge magnetoplasmons (EMP) is presented 
for strong magnetic fields, low temperatures, and integer filling factors v. It is valid for negligible 
dissipation and lateral confining potentials smooth on the scale of the magnetic length io but suf- 
ficiently steep that the Landau level (LL) flattening can be neglected. LL coupling, screening by 
edge states, and nonlocal contributions to the current density are taken into account. In addition 
to the fundamental mode with typical dispersion relation u ~ q x \n{q x ), fundamental modes with 
, acoustic dispersion relation u ~ q x are obtained for v > 2. For v = l,2a dipole mode exists, with 

■ dispersion relation lu ~ q%, that is directly related to nonlocal responses. 

a 

I. INTRODUCTION. 

Edge magnetoplasmons (EMP's) have atracted considerable interest in the past years 0]- |ll[] and Ref. [0 contains 
a discussion of older studies. In JllJ we have treated EMPs, oc A(u>, q x , y) exp[— i(ujt — q x x)], for v = 1,2 and 
very low temperatures when the unperturbed density profile drops sharply at the edge on a length of the order 
of the magnetic length £ . Such a profile is valid for ksT <C hv g /£o, where v g is the group velocity of the edge 
states. Then for v = 1,2 the unperturbed electron density no(y), normalized to the bulk value no, is calculated as 
n Q (y)/n = {1 + $[(y re — 3/)/£q1}/2, where y re is the coordinate of the right edge and <3?(y) the probability integral. 

The main results of Ref. jflj were the existence of many symmetric and antisymmetric modes with respect to the 
edge, some of which can propagate nearly undamped even when the dissipation is strong and were, therefore, termed 
edge helicons. In this quasi-microscopic treatment, however, the contributions to the current density j^iyj, q x , y) were 
obtained from microscopic expressions valid when the components of the electric field E, associated with the wave, are 
smooth on the £q scale. This is not well justified. In Ref. Jll|] we neglected possible nonlocal effects and approximated 
the contributions to j ^(uj , q. x , y) with those obtained when E y (y) is smooth on the scale of Iq. The model jll| also 
neglects the screening of the two-dimensional electron gas (2DEG). Here, using a RPA framework we study the effect 
of nonlocality and of edge-states screening on EMPs for filling factors Q>v>l. For even filling factors v we neglect 
the spin splitting. As we will show, the existence of additional modes for v = 1(2), predicted in Ref. pi] ], is confirmed 
by the present fully microscopic treatment and the results for the fundamental mode remain valid for weak dissipation 

■ despite the neglect of screening and nonlocality in jll|] . 

In Sec. II we present the channel-edge characteristics and derive an integral equation for the inhomogeneous wave 
charge density within a RPA framework. In Sec. Ill we obtain the fundamental and dipole EMP's for v — 1,2 and 
in Sec. IV we consider fundamental EMPs for v = 4 and v = 6. Finally, in Sec. V we discuss briefly the results and 
make concluding remarks. 
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II. BASIC RELATIONS 
A. Channel-edge characteristics 



We consider a zero-thickness 2DEG of width W and of length L x — L in the presence of a strong magnetic field B 
along the z axis. We take the confining potential flat in the interior of the 2DEG (V y — 0) and parabolic at its edges, 
V y = m*£l 2 {y — y r ) 2 /2, y > y r . V y is assumed smooth on the scale of £q — (ft/m*^) 1 / 2 such that SI <C lo c , where 
oj c = \e\B/m* is the cyclotron frequency (e < 0). In the Landau gauge for the vector potential A = (—By, 0,0) the 
one-electron Hamiltonian h is given by 
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h° = [ipx + eBy) 2 + f y ]/2m* + V y 



(1) 



where p is the momentum operator. The eigenvalues and eigenfunctions corresponding to Eq. (|l|) near the right edge 
of the channel, with y$ = yo(k x ) = i"^k x > y r , are well approximated by 

E a = E n (k x ) = (n + l/2)Rw c + m*O 2 (y - t/r) 2 / 2 . (2) 

and 

|a >= ^a(r) = e ik * x V n (y - y Q )/VZ, (3) 

respectively. Here r = {x,y}, a = {n,k x }, ^ n (y) is a harmonic oscillator function. The energy spectrum (|^) of the 
n-th LL leads to the group velocity of the edge states v gn = dE n (k r + ki )/hdk x = hQ 2 ke /m*uj 2 with characteristic 
wave vector k^ = (lu c /HQ) y/2m* Ap n , Ap n = Ep — (n + \/2)Tiuj c , and Ep is the Fermi energy. The edge of the n-th 
LL is denoted by y rn — y r + fgfee = igkm, where k rn = k r + ki , and W — 2y r Q. We can also write v gn = E en /B, 
where E en = Q^/2m* Ap n /\e\ is the electric field associated with the confining potential V y at y rn . We also introduce 
the wave vector k r = y r /(^. 

For definiteness, we take the background dielectric constant e to be spatially homogeneous. Assuming 3> 1, we 

can consider an EMP along the right edge of the channel of the form A(u, q x , y) exp[— i(u)t — q x x)] totally independent 
of the left edge. For v even the spin-splitting is neglected. 



B. Wave charge density and electric potential at the channel edge 

As in Refs. and jL2j, we assume that without interaction the one-electron density matrix is diagonal, i.e., 
< a\p(°'\0 >= f a 5 a p, where f a = 1/[1 + exp((E a — E F )/kBT)] is the Fermi-Dirac function. 

For the application of the RPA we follow the self-consistent field approach of Ref. and references cited therein. 
The one-electron Hamiltonian in the presence of a self-consistent wave potential V(x, y, t) = V(u>o, q x , y) exp[— i(u>ot — 
q x x)j + c.c, is H(t) = h + V(x,y, t). The corresponding equation of motion for the one-electron density matrix p 
reads 

lh f t =[H(t) lP ]-%-p(% (4) 

where [, ] denotes the commutator. On the right hand side (RHS) of Eq. (Q) we have introduced the phenomenological 
infinitesimal term oc 1/r (t — > oo) that leads to the correct rules for contour integration around the singularities, 
cf. Refs. and jl4|]. Notice that r — > oo corresponds to the collisionless case while using a finite r provides the 
possibility of estimating roughly the influence of collisions. Here the most effective collisions are related to intra-LL 
and intra-edge transitions, cf. Refs. [0- (Hi. 

We take the Laplace trasform, with respect to the time t, of Eq. (Q), write R(uj) = / °° e luJt pdt and R a p(uS) =< 
a\R(uj)\(3 >, and look for its solution in a power series in V 



(5) 

n=0 



where R^p(u>) — {if a /u)5 a p. Because we consider linear EMPs, it is sufficient to take into account only the 
first two terms n = and n — 1 on the RHS of Eq. (|J). Then in V(x,y,t) we can consider only the term 
V(^o, q x , y) exp[— i(uiot — q x x)], which leads to 



ff (i), \ _ Kfj - U) < a\V(vo, Qx, yy q * x \(3 > 
a P { ' {w-w Q )[E p - E a + ticj + ih/T] ' 



Taking the trace of p with the electron charge density operator, eS(r — f), gives the wave charge density as 

-oo+i V —p 



/OC' + ?7/ 
dwe-^^iZ^H^W^fr), (7) 
-oo+in _.o 
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where n > 0. From Eqs. (^|) and §j§ it follows that p(t,x,y) — p(t,q x ,y) enp(iq x x). Moreover, for t/r 1, the con- 
tributions related to transitional processes are already negligible. It follows that p(t, q x ,y) = p(^o, y) exp(— iujQt). 
Further, from Poisson's equation the wave electric potential 4>(t, q x , y) induced by p(t, q x , y) is given as 

&(t,qx,y) = - dy'K (\q x \\y-y'\)p(t,q x ,y'), (8) 

^ J — oo 

where Kq(x) is the modified Bessel function; (f> and p pertain to the 2D plane. 



C. Integral equation for EMP's 

For t/r 3> 1 the relation p(t,q x ,y') — p(u>o,q x ,y')exp(—iu>ot) holds; then from Eq. (||) it follows that <p(t,q x ,y) = 
<fr(u)Q, q x , y)exp{—iu>ot). In the absence of an external potential we have V(u)q, q x ,y) = e(j)(ujQ, q x , y). As a result, from 
Eqs. for t/r >■ 1 , we obtain the integral equation for p(u>,q x ,y) as 

p(">^y) = 7L £ E e - £+£ + ih/r n -<^ kxa> M 

n a ,n^— k xa 



X 



oo 



— oo 



dy dy' U, lani3 (y,k xa ,k x g)K (\q x \\y ~ y'\) p(uj,q x ,y'), (9) 



ifi{y,k xo ,,k x p) = ^n a (y-yo{k X a))^n (y-yo(k x g)), k x g = k xa - q x - We dropped the subscript from u . 
For flat LL's, i.e., for O — > and fixed width of the 2DEG, we first apply the Fourier transformation along y to Eq. 
(^) and then carry out the sum over k xa as well as the integral over y' . This leads to p(u>, q x , qy)ei(u>, q x , q v )/e = 0, 
where ti{uj,q x ,q y ) = e;(oj,q) = ei(w,q) is the RPA longitudinal dielectric function (q = {q x ,q y }). That is, for 
p(u>, q x , q y ) ^ we have 

^ = i-%Y. F < > I 2 = o. (io) 

e eqL A ^—J Eg — E a + Twj + in r 

This equation gives the RPA dispersion relation for "bulk" longitudinal wave excitations of a 2DEG in a strong 
magnetic field. It also shows that the "bulk" excitations should have u> > uj c . For definitness, we take u> > 0. 

The integral equation (^) can be considered as a generalization of Eq. (10) of Ref. [|ll| since it takes into account 
nonlocal contributions to the current density oc J dy' a^ ll {y,y')E~ l {y') and the screening by the edge and bulk states 
of the 2DEG. As shown in Ref. [ p!B| , the screening by edge states can be strong. We will obtain only solutions of Eq. 
(|) that are low in frequency, i.e., solutions for which lu/lo c <C 1. The appearence of any wave branch with uj < lo c is 
principally related to the existence of the edge of the 2DEG. Hence, any wave with u> < ui c should be localized at the 
edge and can be called edge magnetoplasmon. 

We consider very low temperatures T satisfying %v gn 3> ioksT. Further, we will assume the long-wavelength limit 
qx^o <C 1, which is well satisfied, e.g., for fundamental EMP |ll]] in the low- frequency region. Then if we compare 
the terms oc fg*, for given ng*, on the RHS of Eq. (0), we obtain that the contribution to the sum over n a with 
n a = ng* is much larger than any other term of this sum or the sum of all terms with n a ^ ng- . The small parameter 
is \lu q x Vgnp* 

(k x g)\/uj c <§; 1, where v gnf1 , (k x g) = h dE n0 , (k x g)/dk x g is the group velocity of any occupied state 
{ng*, k x g} of the ng* LL. We assume that each occupied n-th LL has one intersection with the Fermi level at the edge 
of the channel and denote the group velocity of its edge states as v gn = v gn (k rn ). The small parameter given above 
implies q x v g o/uj c <C 1, as v g Q is typically the largest among v gn . Similar results follow from an analysis of the terms 
oc f a * m the sum over ng on the RHS of Eq. ([)]). Hence, for w <C oj c and q x v g o <C ui c the terms with n a ng can be 
neglected. This leads to the integral equation 

p(u, q x , y) = ^— V V — ^"'^"p' f n «< k *<* i ■%/ n «««= (f ' fc ^> k ™ _ Q*) 

poo pOO 

x dy dy 1 II 

- qx)K (\q x \\y - y'\) p(w,q x ,y'), (11) 



where n denotes the highest occupied LL. For even v the RHS of Eq. (|llj) should be multiplied by 2, the spin 
degeneracy factor. We now study EMPs following from Eq. (O). 
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III. EMPS FOR v = 1(2) 



We first consider the case v = 1 and then explain how the pertinent formulas should be modified for v = 2. For 
v = 1 we have n = and Eq. ( pi]) takes the form 

p(ui,q x ,y) = —— \ dk xa °J — _2_*« n 00 (2/, t„, fc xa - q x ) 
■KheJ_ 00 u-v g o\k xa )q x 

x dy dy' n o(y, \\y - 2/1) p(u,q x ,y'), (12) 



where Co = u> + i/r. As in Ref. 11 1, we seek a solution of Eq. ( pL 2|) in the form 

oo oo 

po(u, q x ,y) = XI Po^K Qx)H n (y/e ) = p n \oj, q*)*n{v)Mv)- (13) 

n=0 ra=0 

Here y — y — y r o = y — yoikro) and Eq. ( |l3| ) is the ea;ac£ solution due to the orthogonality of the Hermite polynomials 
H n (x). The expansion fll3| ) is specific to v = 1(2) under the assumed conditions. This implies that the charge 
distortion far from the LL edge, > 10£o, is much smaller than near this LL edge, < to. Physically this assumption 
can be well justified. We call the terms n — 0, 1, 2, the monopole, dipole, quadrupole, etc. terms in the expansion 
of po(co,q x ,y) relative to y = y r Q. Notice that expansion ( |l3| ) can also be understood as an expansion in a complete 
set of oscillatory wave functions corresponding to all, occupied (n = 0) and empty (n > 1) LLs. 

We now multiply Eq. ( [12] ) by H m (y/£o) and integrate over y from y r Q — Aj/o to y r Q + Ayg, where Ayo > 3£q. With 

the abbreviation p^ l \uj, q x ) = Pg™ 1 ', we obtain 

=^L dk ~ ^-Z^°£: *-(*.«-)g(S) 1/,, ffi(*,i*i.«-)pS fl , (»> 

where 5k xa = k xa - k r0 , 5k xa ' = 5k xa - q x , 

1 f°° 

d m0 (qx,Sk xa ) = —== dy H m (y/£ ) U o(y,Sk xa ,Sk xa ') (15) 
V2 m m! J-oo 

and 

b { V(q x ,\q x \,6k xa )= dy dy' U 00 (y, 5k xa , 8k xa ,)K (\q x \\y - y'\)* n (y')* (v') ■ (16) 



In both d m o(q x , 5k xa ) and b^ \q x , \q x \, 5k xa ) the first argument represents the ^-dependent term of the argument of 
the wave function "$>v{y — 5k > ) and, if replaced by 0, it means this term can be neglected. We will also use the 
coefficients a mn (q x ) given by jM 



a mn (qx) = a nm (q x ) = dx ^ m (x)^ (x) / dx' K (\q x \\x - x'\) ^ n {x')^ {x'), (17) 
and satisfying &£o(0, \q x \,ti) = a n0 (q x ). 



A. Fundamental EMP 



Let us approximate the numerator on the RHS of Eq. (|TJ) by the first nonzero term in the expansion over q x : 
fo,k xa -q* ~ /o,fc XQ ~ -qx{df Q ^ a /dk xa ) = S(k xa - k rQ )q x - Then Eq. (|lj), after integration over k xa , gives 

P { o l) = 4- e | d mO (lx,0)± (|^) 1/2 eVl^O) &\ (18) 

n=0 
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where ljq = (to — v g oq x ). Further, if we neglect a small term oc q x £o in the argument of the wave function ^>o(y + £qQ x ), 
we obtain d m0 (q x , 0) w <i m0 (0, 0) = S m . and &„ |g»|,0) « &i° (0, |?x|,0) = a„ (fe)- Substituting the former in Eq. 
(|l8|), we obtain p^ = for m > 1; thus, only Pq "* can be finite. Then Eq. (jl^) can be rewritten as 

C°) e2 qx i \ (°) do^ 
Po = —r — aoo{qx) p '■ (19) 

irne LUq 

For Pq -* 7^ Eq. (|l^) gives the dispersion relation 

2 1 3 i 

w = q x v g a + -<j° q x [\n(- — — ) + -] - -, (20) 
e |<Zx|«0 4 r 

where = e 2 /2itTi for = 1 and cr° = e 2 /nh for z/ = 2. We have used the result a o(<fe) « [ln(l/|g x |^o) + 3/4] 0] 
for \q x \^o *C 1. If we neglect the dissipation, the dispersion relation ( pp| ) and the corresponding charge density profile 
coincide with the corresponding expressions for the fundamental mode obtained in Rcf. jyj. 

Untill now nonlocal effects and screening by the 2DEG have not changed the results for the fundamental EMP of 
Ref. . In Ref. (llj it is shown that the next term that can affect the fundamental mode of the n = LL is the 
quadrupole term. Then, if we neglect dissipation ui — q x v g o becomes slightly larger, by a factor [l + 0.125/<2oo(<Zx)] ~ !• 
Typically we have aoo(q x ) > 10. In addition, the spatial structure of the fundamental EMP acquires a small quadrupole 

term \p^ /p^\ ~ l/8a Q0 (q x ) <C 1 JXl|] ■ It can be shown from Eqs. (p^)-(p^) that here too the corrections to the 
fundamental EMP due to the quadrupole term are very small; namely, aoo{q x ) on the RHS of Eq. (|o|) should be 
changed to 0,00(92;) + G2o(Sx)9x^o/2"v2 ~ aoo(qx) — q^o/S. Further, the spatial structure of the fundamental EMP 
acquires a small quadrupole term: p$ / p^ = q x l 2 /& <C 1. Thus, for weak dissipation the results of Rcf. [J^TJ] for the 
fundamental EMP are nearly the same as those of the present microscopic treatment. 



B. Dipole EMP 

For a more acurate calculation of the RHS of Eq. ( |l4| ) we must take into account higher order terms in the expansion 
over q x . As shown below, when this is done it leads to additional branches. To obtain the dipole branch correctly we 
must keep the first three nonzero terms in the expansion of the numerator of RHS of Eq. ( J14| ) over q x by writing 

fo,k Ka -q x ~ fo,k xa ~ [Qx ~ h ^g^2~^ S ( kxa ~ kr °)- ( 21 ) 

Further, we will consider only the first two terms, n — and n = 1 in the sum of Eq. (|l^) . Then from Eq. ( |l4| ) for 
in = and m = 1 we obtain the following system of linear equations for p °^ and p^ 



(0) e 2 
Po 



nheuiQ J_ 



1 /2 

dk xa (fo.k 3:a -q x - fo,k xa ) d 00 (q x ,5k xa ) ^ (2 n nl\ b^ \q \,Sk xa )p ( n \ (22) 



n=0 



lJ) -— I dk xo ,(f 0jklca - qx - f ,k xa ) d w (q x ,Sk xa ) ^ (|rjj) b ( ° \q x ,\q x \,5k xa )p n) . (23) 



irheCbo 



We assumed that, for T — > 0, v g o(k xa ) is independent of k xa in the vicinity of k xa = k r Q and replaced it by v g o- This 
assumption does not impose any important restriction on the analysis. Therefore, taking into account Eq. (Ell) we 



can replace 1/ ((I) — v g o(k xa )q x ) by l/(w — q x v g o) and place it in front of the integrals in Eqs. (22)-(E3|). 

In Eq. (G3) we can assume doo(q x , 6k xa ) = 1 since corrections to it, varying as ~ (q x £o) 2 , are not essential. To 

evaluate the contribution oc p °^ on the RHS of Eq. J22]) it is sufficient to keep only the first term, oc q x , on the RHS 



of Eq. (J21J). This gives 

„2 



Po Q) = ^^feC (0, l&|,0)4 0) + V / 2J ( 1 ) (gx)p^ ) }, (24) 



where 
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I oi(Qx)= dkx a (fo,k xa -q x - fo,k xa )bio(q x ,\Qx\,Skxa), (25) 

J — oo 

and where b \q x , |<?x|>0) has justifiably been approximated by 6^(0, |<7x|>0). Before calculating ioiH<7x) we notice 
that the approximation for doo(q x ,6k xa ) entails the neglect of small corrections to I^{q x ) proportional to (q x £o) k if 
k > 4. Then Eqs. @ and Eq. @ give 

^oPfe) = y^-fe4) 2 aii(9x) + 

q 2 f°° q d 9 to 

+lf / dk xa 5{k xa - fc r0 )[l + ^^7— ]^T— & io (fej |fe|,^xa)- (26) 

* J —oc " OK xa OK xa 

The corrections oc (q x £o) k > k > A, we neglected in the first term on the RHS are related to the expansion 

*o(y + ilq x ) = MV) - + \[V2^ 2 (y) - * (y)](g^o) 2 + ... (27) 

It can be shown that the term (q x /3)d/dk xa in the square brackets [...] of Eq. (|26| ) gives, upon integration over 
k xa , a vanishing contribution since the remaining double integral changes sign upon making the changes y —y 
and y' — > — y' . As for the term corresponding to 1 in [...], its evaluation gives q x £oau(q x )/V^- Hence, the two finite 
terms on the RHS of Eq. (|26| ) cancel each other. The final result is Iq\ (q x ) — and Eq. ( pi] ) takes the form of Eq. 
©. That is, in this approximation the fundamental mode is purely monopole and is totally independent of dipole 
excitations. 

We now consider Eq. (p3|). Here we should consider six different contributions I^(q x ), (<7 X )) and lj„ (q x ) for 
n = 0, 1 related, respectively, to the first (oc q x ), second (oc q x ), and third (oc q x ) term on the RHS of Eq. (|2l]). Then 
Eq. (Ba) can be rewritten as 



rf'-^EE^vf. (») 



where 



u n=0 (=1 



J&Ofe) = (^J ) &dio(&,0)^(gx,kx|,0), (29) 
■*in (&0 = ^ " 9 2 / dk xa d w (q x ,5k xa )b { °l(q xi \q x \, Sk xa )— — 5(k xa - fc r0 ), (30) 

*• J —oc (J& X u 

^ln (&) = 51 * 9x / ^ Q rflo(<7 K ,^a)&i°o ) (feJ'? K |,^fcK Ce )^73- (5 ( fc 2;a ~ k ro)- (31) 

Again we will neglect contributions oc (q x £o) m for m > 4. In Eq. ( p9| ) we have diofe, 0) = — (q x £o) / V^+O((q x £o) 3 ) 
and 6^(^,19,1,0) = a ao (q x ) + O((q x £ ) 2 ), b^(q x , \q x \, 0) = -(&^ /V2) ttu( fc ) + 0(( fe 4) 3 ). It follows 

^io (fe) = 7Tn — a 00 (q x ), In(q x ) = — ^ — a n (q x )- (32) 

Evaluating the integral in Eq. (|3^) by parts we obtain 

= (2 "' 1 2 f )V2 ql [dio(q x ,6k xa )] lskxa = b^(q x ,\q x \,0) + 

dw(q x ,0) -^j— [ & io(9x, \q x \,Sk xa )]\5k* a =oY ^ 

This gives the contributions 
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Uo (?x) = — aoofe), ^li (<?*) = 7 anfeO- (34) 



Evaluating the integral in Eq. (plh by parts we obtain 



/<?(«.) = 0, ig>( fc ) = oufe). (35) 



Using Eqs. (|32|) , ( 34 ) and (|35|) we can rewrite Eq. ([281) as 



(i) _ e^oufe) 3 (1) 

Eq. ( |36| ) shows that the dipole branch, similar to that in Ref. , is not coupled to monopole excitations of the charge 
density. Further, Eq. (|36|), with an(q x ) — 1/2 for \q x \£o <§C 1 and precision < 0.2% pJ| , gives the dispersion relation 

of the dipole branch, Pq (u,q x ) 0; 

u = v g0 q x + —a°ql , (37) 

oe " t 

where az x = e 2 /2irh for ^ = 1 and a ~ = e 2 /ttH for ^ = 2. 

An important difference of Eq. ( p7| ) from the corresponding result, Eq. (40), of Ref. for the pure dipole EMP 
is that here we have u) — v g oq x oc q x whereas in Ref. [|l If] this becomes u> — v g oq x cx g^. That is, the term caused by the 
Coulomb interaction has a nonacoustic behavior if we neglect dissipation. The calculations demonstrate clearly that 
the spatial dispersion effects, although they depend on q Xl they are essentially related to the wave structure along the 
y-direction. It is known that spatial dispersion is directly related to the nonlocality of responses, as expressed by the 
dielectric function, et c. [jl6| . In Eq. (j37 ) both nonlocal effects and edge states screening are taken into account. We 
point out that in Ref. ]17| a similar dispersion law, u> cx aq x , was obtained for short-range forces, i.e., for an interaction 
that is essentially not of the Coulomb type. 

IV. FUNDAMENTAL EMPS FOR v = 4 AND 6 

For v — 2(fi + 1) the n = 0,...,rt = n LLs have intersections with the Fermi level at y r o,..., y r n, respectively, and Eq. 
( |il"| ) can be written as 

p(u,q x ,y)= y~" [ dk x J n ^ a ~ q * ^"- fc *° U nn (y, k xa , k xa - q x ) 
€ Z~oJ-°° u — Vgn{k xa )q x 

/OO />00 
dy j dy' U nn (y,k xa ,k xa - q x )K (\q x \\y- y'\)]p(uj,q x ,y'). (38) 
-00 J — 00 

For confining potentials smooth on the £q scale we have Ay m _i jm = j/ rm _i — y rm 3> £q, where m < h. For v — 4 there 
is only one inter-LL length Aj/o.i = y r o ~ Uri- Further, making the same approximations, in the long- wavelength limit 
q x £o "C 1, as in Sec. IIIA and integrating over k xa in Eq. (|38|), we obtain 

9 2 ™ 

ze T 2 , 



P(u,q x ,y) = —r- Y) ^-^l(V - Vrn) 



n=0 



x dy dy'[^ 2 n (y-y rn )K (\q x \\y-y'\)]p{u;,q x ,y'), (39) 



where u) n = u)—v gn q x . It follows that p(u>, q x ,y) can be represented by a sum of charges p n (u>, q x , y) — p„ (u), fe)^ 2 {y— 
y rn ), localized at the edge of the n-th LL, within a region of extent ~ y / 2n + 1£q, n = 0, ...,n. The result is 

n 

p{u;,q x ,y) =^p!£\u),q x )*l(y-y rn ). (40) 

n=0 

We substitute Eq. ( |40| ) into Eq. ( |39"|) and demand that the coefficients of 'J 2 (y — y rn ) on both sides of Eq. ( |39| ) be 
equal. This leads to fi + 1 linear homogeneous equations for p^\^ , q x ) ■ 
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A. v = 4 



Due to Ayo.i 3> £o we can neglect the exponentially small overlap between po(u>,q x ,y) and p\(uj,q x ,y). We will 
assume g x Aj/oi 1. As shown above, from Eqs. ( j39| ) and (f4(i|), with h = 1, we obtain the system 

= ^^ [aoofe)p (°) + all(q x ,Ay 01 )pM}, (41) 
(0) = 2^fe n ( 0) +a H fe)0 ) p (°)], (42) 

where pn\w,q x ) = pn and 

/>OC />oc 

a™™(q x ,Ay) = / cfe (i)*J,(i')ifo(|g»| \x - x' + Ay\). (43) 



J — oo */ — oo 

Here a%™(q x ,Ay) = a™(q x ,Ay\ a™™(q x ,Ay) = a™"^, -Ay) and agg(fe,0) = a 00 (q x ). For a nontrivial solution 
of the system of Eqs. (Eh) and M2) the determinant of the coefficients must vanish. This gives the dispersion relation 
of the renormalized fundamental EMP of the n — LL as 

u = q x voi + ^ fePMirV) - M^ 1 ) + |] - (44) 
e IfeKo «o 5 t 

and that of the n = 1 LL as 

. n ,Ay ix . 2 i 

w = fevoi + — ?x ; — ) + v ) (45) 

e to 5 r 

where cr° x = 2e 2 / tt?z. and Wqi = (v g o + v g x)/2. All double integrals involved in the coefficients of Eqs. ( [4l| ) and 
(E2n can be evaluated within the approximation Kq(x) rs ln(2/|x|) — 7, where 7 is the Euler constant. We obtain 
oJ5(?x,0) « ln(l/|g x |4) + 1/4 and ajjfe, Ay i) « (m(2) - 7) - ln(| fe | Ay i) ss ln(l/|g x Ay i|) + 0.1. Notice that in 
the absence of dissipation Eqs. ( p4| ) and ( p5| ) coincide with the quasi- microscopic results of Ref. [Q. ^,From Eqs. 
( |4l| ) and ( |42| ) it is easy to see that, if the inter-LL Coulomb coupling is neglected by setting aH(q x , Aj/oi) = 0, the 
dispersion laws of the decoupled fundamental EMPs of the n — and n = 1 LLs are given, respectively, by Eq. (pf~ 
and by 



2e 2 , 1 , 1. i 
-Tlx Mt-jt + I - -■ 



w = feWgi + — — fe[hi(— —) + -] - -. (46) 



Substituting Eq. (gj) in Eq. @ and Eq. (|f) in Eq. @ we obtain, respectively, p^/p^ » 1 and p^ / pf> » -1. 
This means that in the former case the wave charges localized at the edges of the n = LL and n = 1 LLs are in 
phase whereas in the latter they are out of phase. Therefore, the former EMP has an acoustic spatial structure along 
the y axis while the latter EMP has an optical spatial structure though its dispersion law is purely acoustic. 

B. v = 6 

For v = 6 we have the intersection of the n = 2 LL with the Fermi level, at y r 2, in addition those of the n = and 



n = 1 LLs; n = 2. Corresponding to Eqs. (39) and Eq. (Wfj) we now obtain 



Po 0) = ^^[aoo(q x )Po 0) +aoo(9x»At/oi)Pi 0) + aoofe, At/ 2)p 2 0) ]> (47) 
7rn,e wo 

pi 0) = ^Iqa^Ayoi)^ +a\l(q x ,0)pf ) + aH(q x , Ay 12 )$\ (48) 

pi o } = 2e^|x [a 22 fe)A2/o2)p W +a??fe;Ayi2)p (o) +^(^,0)^°)]. (49) 

7Tfte W2 
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The vanishing of the determinant of the coefficients leads to the cubic equation 



to 



'3 



+ a 2 {q x )u>' 2 + a\{q x )w + a Q (q x ) = 0, (50) 



where u> — uj / (2e 2 q x / nhe) . The expressions for the coefficients a,k(q x ), k = 0, 1, 2, are given in the appendix. We will 
assume that q x Ayij <C 1, i ^ j < 2. 

It can be shown that all three roots of Eq. (|o|) are real and different as they correspond to the irreducible case. 
Assuming Ay 2, Ayi 2 , Ay i 3> ^o, these roots give the dispersion law for the renormalized fundamental EMP of the 
n = LL as 

, 2a yx row 1 ^ 2 n i^V^\ , i Avm\ , , / A J/12m , 8 i i 

u = feVo12 + ^r^ {31n( ^ } - 3 [ln( ^r } + ln( ~ } + ln( ~ )] + is } - ^ (51) 

and that for the renormalized fundamental EMP of the n = 1 LL, w + , and n = 2 LL, w_, as 

. 2^ n , Ay 02 Ay i ^ Ay 12 7 

«± = <» 2 + _ fc {]n(-^-) + ln( — ) + ln(_) + -} ± 

^^{ln 2 (^) + ln 2 (^) + ln 2 (^i) + A}V» - I. (52) 
9e Ayoi Aj/12 Ayi2 r 

Here = 3e 2 /nh, v i 2 = (v g o + v g \ + v g2 )/3, and 

A = (v g0 + 3/4) [(v g0 - v g2 + 3/4)/2 + ln(Ay 02 Ayoi/Ay 12 2 )} 
+{v g i + 1/4) [(^i - v gQ - l/2)/2 + ln(Ay iAy 12 /A 2 / 02 2 )] 

+v g2 [(v g2 - v g i - l/4)/2 + ln{Ay 02 Ay 12 /Ay 01 2 )}, (53) 

where u 9 i = v g i/ (2e 2 /nhe), i — 0,1,2. Notice that the fundamental EMPs of the n = 1 and n = 2 LLs have 
purely acoustic dispersion laws, cf. Eq. (|5^) . Only the fundamental EMP of the n — LL, Eq. (|Tj), behaves in the 
previously obtained manner oc q x ln(q x ), cf. Refs. jO- JhJ), for a fundamental EMP. It is difficult to clearly disentangle 
the contributions of the n — 2 LL in Eq. (p2]). It can be shown though that, if they are neglected by setting a 22 , a 2 f, 
fl oo; v g2 to zero, then we have w 2 — and Rec<j_ = as we should in this limit. 

The solutions of the cubic equation, Eqs. ( |5l| ) and (|52|), are not obtained using the standard expressions, reproduced 
in the appendix, which become very unwieldy for the present case. Instead, the following reasoning is used. If uj has 
an acoustic character for some branches, then uj in Eq. (^0|) is independent of q x and should satisfy the quadratic 
equation obtained by differentiating Eq. ( |50| ) with respect to q x 

_g(gx) < 2 . dai{q x ) , da (q x ) _ 

— ^ OJ H q Ul H — U. (04J 

oq x oq x dq x 

The two roots u) 2 and W3 of Eq. (|w|) are indeed independent of q x and lead to the solutions lj_ and cj + , respectively, 
given by Eq. (|52|). The other root is given by uj\ — —a 2 (q x ) — (oj 2 + lu 3 ) and results in Eq. (|5l]). Of course both 
methods give the same solutions w,-, j = 1,2,3. Notice that on the RHS of Eqs. ( [5l] ) and ( |52] ) we used approximate 
expressions of a J J(q Xl Ayij) which , assuming Ay^jl^ 3> 1, i ^ j, are obtained with the approximation 

alx dx'*1(x)* 2 {x')\n(\x - x' + A VlJ \/l ) a ln(A^/£ ). (55) 



— 00 j — 00 



Below we will demonstrate that the approximation ( ]55j ) is well justified even for not-too-large values Ay^jl^. For 
that we will calculate the dispersion laws of the renormalized fundamental EMPs using the exact expressions of 
a{f (q x , Ay^), i.e., with the double integral on the LHS of Eq. (|55| ) evaluated numerically. In the latter case we will 
calculate ujj from the standard expressions (|5^) -(^lj) of the appendix. 

For a GaAs-based 2DEG and negligible dissipation the dispersion laws (|5l|) and ( |52| ) for the renormalized, by 
the inter-LL Coulomb coupling, fundamental EMPs are shown in Fig. 1 by the solid curves. The parameters are 
m* w 6.1 x 10~ 29 5, e w 12.5, and S]«7.8x 10 11 sec _1 H). For f = 6 and B = 3Tesla these parameters lead to 
%uj c w 5meV, lu c /Q sa 10. Here = 2e 2 /irheio is a characteristic frequency almost equal, for the conditions stated, to 
lu c . We have also assumed Ap 2 — hui c /2. This gives v g o w v5^0j u gi ~ v3fi^0i w g2 ~ ^^0, Aj/02/^0 ~ (w c /f))[\/5— 1], 







~ (<^c/^)[v / 5 — VS\, and Ayyil^-a ~ (w c /£2) [V3 — 1]. The dashed curves show the fundamental EMPs of the 
totally decoupled n = (Eq. (|2(i|), topmost dashed curve), n = 1, and n = 2 (lowest dashed curve) LLs, i.e., the 
modes obtained by neglecting the inter-LL Coulomb coupling. In this case all coefficients a 3 J on the RHS of Eqs. 
(p7|)-(p9|) vanish for i ^ j. This leads directly to the dispersion laws represented by the dashed curves in Fig. 1. 
The dotted curves show the dispersion laws of the renormalized fundamental EMPs using the exact expressions for 
a" (q x ,Ayij), i ^ j. It is seen that each dotted curve is very close to the corresponding solid curve. 

For the same parameters and conditions as in Fig. 1 we show, in Fig. 2 , p(u> , q x , y) / p* (solid and dotted curves) 
and p(u>,q x ,y)/2p* (dashed curve) as a function of Y = (y — y r 2)/^o with p* — p °\uj,q x )/^/ir£o. The solid curve 

corresponds to the topmost solid curve in Fig. 1, Eq. (|5l|), with Pi°Vpo°' ~ 1-0 and p% / p^ » 1. As can be seen the 
renormalized fundamental EMP of the n = LL has an acoustic spatial structure along the y axis. The dashed curve 
corresponds to the lowest solid curve in Fig. 1, cj_ of Eq. (|52]), with Pi / Pq ~ —2.0 and p% / Pq ~ 1- Hence the 
renormalized fundamental EMP of the n = 2 LL has an optical spatial structure though its dispersion law is purely 
acoustic. The dotted curve corresponds to the middle solid curve in Fig. 1, uj + of Eq. (p2[), with p\ / Pq ~ —0.1 and 
/ 9 2°Vpo°' > ~ ~l-0. Thus, the renormalized fundamental EMP of the n = 1 LL has an optical spatial structure though 
its dispersion law is purely acoustic. Notice that the dependence of p^/Po^i j = 1, 2, on q x is typically weak. 

In Fig. 3 we use the same values for m*, e, fl, and Ap 2 = hoj c /2, v = 6 as in Fig. 1. However, we take B = 1.5Tesla 
and this leads to u> c /£l w 5. This gives y/2 times larger values for v g o » \/5Q£o, v g \ \/3^0: and v g2 « fl£o> and 
twice smaller values for Ay 02 /£o ~ (w c /fi)[v5 — 1], Ay i/£ ~ (a; c /fi)[V5 — V3], and Ayi 2 /£ ~ (w c /fi)[\/3— 1] as 
compared with those in Fig. 1. All curves are marked as in Fig. 1 and show the corresponding dispersion laws. 

In Fig. 4 we plot p{oj , q x , y) / p* as a function of Y = (y — y r 2 ) / f° r the conditions corresponding to the solid curves 
in Fig. 3. All curves are marked as in Fig. 2. We have, respectively, p^/p^ « 1.0 and p 2 °^ j Pq ' ~ 1 (solid curve), 

Pi / Pq ~ — 2.1 and P^Vpo ~ ^-^ (dashed curve), and Pi" 1 / p^ ~ —0.1 and P2 I Pa ~ —1-0 (dotted curve). As in 
Fg. 2, the renormalized fundamental EMP of the n = LL has an acoustic spatial structure along the y axis while 
those of the n = 1 and n — 2 LL have an optical spatial structure though their dispersion laws are purely acoustic. 
As noted above, the dependence of Pj /Po°\ j — lj 2, on q x is typically weak. 



V. DISCUSSION AND CONCLUDING REMARKS 



We presented a fully microscopic model of EMPs in a RPA framework valid for integer v > 1 and confining potentials 
that are smooth on the £q scale but still sufficiently steep that LL flattening can be neglected 20 . The model takes 
into account LL coupling and treats only very weak dissipation. The main results of the present work are as follows. 

i) For moderately steep confining potential we presented a microscopic model that improves the quasi-microscopic 
approach of Refs. (ll]] and flisfl . In particular, the model combines features of the following distinct edge wave 
mechanisms. In mostly classical models, e.g., pi, [Q, the position of the edge does not vary but the charge density 
profile at the edge does. In a sence this edge wave mechanism is the analog of that for the Kelvin wave |H[ at the edge 
of a rotating "shallow" sea with chirality determined by the Coriolis parameter which corresponds to the cyclotron 
frequency oj c . Another edge-wave mechanism, fully quantum mechanical, is that of Refs. p~(i|| , |^7| , and (||]- p3[ , in 
which, for v = 1, only the edge position of an incompressible 2DEG of the lowest LL varies and with respect to which 
the density profile is that of the undisturbed 2DEG. The approach of these works is limited to the subspace of the 
lowest LL wave functions, neglects LL mixing and dissipation, and results in a single EMP with dispersion law similar 
to that in Ref. §. 

ii) We confirmed, for v = 1(2), the existence of EMP modes in addition to the fundamental EMP and single one 
obtained for vertically steep unperturbed electron density profile 0] . This is in line with our earlier quasi- microscopic 
results |ll[] for moderately steep [jllj confining potentials. The additional modes result from an exact solution of Eq. 
( p"3| ) in terms of the complete set of the Hermite polynomials. With this expansion we can also make contact, for 
v = 1(2), with other microscopic theories JlCj] , fTi| , and |p2fl - ]23| that are limited to the subspace of the lowest LL 
wave functions. If we retain only the n = term in Eq. (]13[), we obtain only the n = LL fundamental mode with 
the same dispersion relation oj ~ q x h\q x . The additional modes result from retaining ( correspond to) higher-order 
terms in the expansion Eq. (|l3|). This shows that the limitation to this LL subspace is too restrictive and indicates 
the importance of LL coupling even for v = 1(2). 

iii) An important additional mode is the dipole EMP, presented in Sec. Ill B, with dispersion relation uj ~ q x . This 
differs markedly from the corresponding result lu ~ q x of Ref. JlT[ and signals nonlocal responses that were previously 
neglected. This and other differences between the present results and those of Ref. [[flj further indicate the importance 
of having a fully microscopic theory for EMPs. 
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iv) As we showed, taking LL coupling into account is an essential ingredient especially for v > 2. As Fig. 1 and 
Fig. 3 demonstrate, the coupling strongly renormalizes the uncoupled fundamental LL modes and results in one mode 
behaving as the usual one uj ~ q x In q x and the others behaving in an acoustic manner u> ~ q x . This partly reminds the 
results of Ref. but the resemblance should not be overestimated because the two models are drastically different. 

v) We have not treated spin effects beyond the HA or RPA. Our study is focused on important wave effects of 
non-spin nature and spin splitting is assumed negligible for even v. Thus, skyrmions and spin textures are not dealt 
with in this work. Though neglecting the spin splitting is a reasonable approximation, for even v > 2 in the bulk of 
the channel, its validity near the edges remains uncertain in view of the work of Refs. [p"5| and 

(vi) Finally, a few remarks are in order about the studies of Refs. |2a] and ggj that used the RPA. The study of 
Ref. |^6| is too simplified and in essence repeats the results of Ref. |25| and of works cited therein for the " optical" 
EMP. As for the results of Ref. pq], due to omitted important contributions, they are essentially different from ours 
for both the acoustic and the optical modes. As Sec. IV A shows, an important logarithmic term, caused by the 
Coulomb interaction, is missed in Eqs. (54a) and (54b) of Ref. j25j. In addition, Refs. Q and did not study a 
dipole mode or multipole modes pertaining to any occupied LL. These modes of ours are totally different than those 
of Sees. 3.1.6 and 3.1.7 of Ref. M. 
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APPENDIX 

The expressions for the coefficients ak(q x ), k = 1, 2, 3 in Eq 



J) are as follows: 
a 2 (q x ) = -[Soo(&) +aii(?x,0) 0)], 



(56) 



and 



a%(q x ) = aoo(fe)aii(fe,0) + a Q0 (q x )aH(q x ,Q) + a\\(q x , Q)a 2 2 \(q x , 0) - 
[<${q x , Ay 02 )] 2 - [aJjGfe, Ayoi)] 2 - [a 22 ^, Ay 12 )} 2 , 

a (q x ) = a o(fe)[a 2 i(fe, Ayi 2 )] 2 + a{{(q x , 0)[aH(q x , Ay 02 )} 2 + aH(q x , 0)[aH(q x , Ay i)f 



(57) 



aoo(qx)a 11 (q x ,0)d 22 (q x ,0) - 2a 00 (q x , Ay i)a n (q x , Ay 12 )a 00 (q x , Ay 02 



(58) 



Here a o(qx) = a o(q x ) + v g0 /(2e 2 /nhe), and a\l(q x ,0) = al\(q x ,0) + v gl /(2e 2 /Trhe), i = 1,2. We also have 
aM(q x ,0) ~ \n(l/\q x \£ ) and, for Ay 02 ,Ay 12 > £ , a 22 (q x ,Ay 02 ) « \n(l/\q x \£ Q ) - \n(Ay 02 /£ ) + 0.1, and 
a 22 {q x , Ay 12 ) » ln(l/|g x |4) - ln(Ay 12 /£ ) + 0.1. 

Eq. (p0|) corresponds to the irreducible case Q{g x ) = q 3 +r 2 < 0, where q — [a\ — a|/3]/3, r = (aia2 — 3ao)/6— 02/27. 
With p — y/r 2 — Q and 8 = arctan(^/— Q/r) the roots obey + uj 2 + — —a 2 and are given by 



uo\ = 2p 1/3 cos(6»/3) - a 2 /3, 



(59) 



and 



uj 2 = -p 1/3 cos(0/3) - a 2 /3 - V3p 1/3 sin(6'/3), 



lo 3 = -p 1/3 cos(0/3) - a 2 /3 + V3p 1/3 sin(6»/3). 



(60) 



(61) 
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FIG. 1. Dispersion relation for v — 6 and B = 3Tesla. The dashed curves are, from top to bottom, the decoupled 
fundamental modes of the n = 0, n = 1, and n = 2 LLs, respectively. The solid and dotted curves are for correspondingly the 
same coupled modes, uj, — 2e 2 /irhedo. For the solid curves approximation (Eq. ( p5|)) is used. 

FIG. 2. Spatial structure of the renormalized fundamental modes of the n = 0[ solid curve), n = 1 (dotted curve), and 
n — 2 (dashed curve) LLs, whose dispersion laws are shown by the solid curves in Fig. 1. The solid and dotted curves show 
p(cu, q x , y)/ p* and the dashed curve p(uj, q x ,y)/2p* as a function of Y = (y — y r 2)/£o with p* = p^\lo, q x )/\^io- 
FIG. 3. Dispersion relation for v = 6 and B — 1.5Tesla. All curves are marked as in Fig. 1. 

FIG. 4. Spatial structure, p(u),q x ,y)/p*, of the renormalized fundamental modes of the n = (solid curve), n = 1 (dotted 
curve), and n = 2 (dashed curve) LLs, whose dispersion laws are shown by the solid curves in Fig. 3. 
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Fig. 1, Balev et al. 




Fig. 2, Balev et al. 




Fig. 3, Balev et al. 
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